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Abstract 



The effect of density inhomogeneity on the YORP effect for a given 
shape model is investigated. A density inhomogeneity will cause an 
offset between the center of figure and the center of mass and a re- 
orientation of the principal axes away from those associated with the 
shape alone. Both of these effects can alter the predicted YORP rate of 
change in angular velocity and obliquity. We apply these corrections to 
the Itokawa shape model and find that its YORP angular velocity rate 
is sensitive to offsets between its center of mass and center of figure, 
with a shift on the order of 10 meters being able to change the sign of 
the YORP effect for that asteroid. Given the non-detection of YORP 
for Itokawa as of 2008, this can shed light on the density distribution 
within that body. The theory supports a shift of the asteroid center of 
mass towards Itokawa's neck region, where there is an accumulation of 
finer gravels. Detection of the YORP effect for Itokawa should provide 
some strong constraints on its density distribution. This theory could 
also be applied to asteroids visited by future spacecraft to constrain 
density inhomogeneities. 



1 Introduction 



The YORP effect for asteroids, initially proposed by Rubincam [9], has now 
been verified 0C1UIH] and is being invoked to explain a variety of evolutionary 
dynamics and states for the smaller members of the asteroid population. 
Despite this success, the YORP effect for the one small NEA whose shape 
and mass is best characterized, the asteroid Itokawa [HI BJ [Q, has not been 
detected yet, despite predictions that YORP should have been evident by now 

Hn]Q The non-detection of the YORP effect for Itokawa raises a number of 

1 It should be noted that a detection of YORP for Itokawa was reported in [6], but was 
subsequently retracted by that author at the 2007 DPS meeting. 
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interesting modeling questions about the details of the YORP effect. 

Rubincam noted that density inhomogeneity could affect computed YORP 
torques [§], although this has not been explicitly investigated in the litera- 
ture. In the following we will show that density inhomogeneity within an 
asteroid has a direct influence on the predicted YORP spin rate evolution, 
and present the theory to compute this effect. We present computations of 
the change in the YORP effect for high-precision shapes of Itokawa. We 
find that center of mass shifts on the order of a few meters can significantly 
change the YORP spin-rate deceleration, and that shifts on the order of 10- 
15 meters could change its sign and make Itokawa's spin rate accelerate. The 
current non-detection of YORP for Itokawa is consistent with a density con- 
centration in the neck region of Itokawa. Once the YORP effect for Itokawa 
is detected, the theory presented here can help us constrain its density in- 
homogeneity. Similar YORP measurements for asteroids to be visited in the 
future could also provide insight into their density distributions. 

2 Theory 

At any given instant the force and moment acting on an asteroid due to in- 
cident, reflected and reemitted solar radiation from its surface can be repre- 
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sented as a function of the sun's location in an asteroid-fixed reference frame. 
In [11] the representation of these forces and moments using a Fourier series 
expansion is described in detail, and algorithms for the computation of the 
coefficients, starting from a polyhedral shape model, are given. The direction 
of the sun in the asteroid body-fixed frame can be described by a unit vector, 
u. For a given solar location we can compute the total force and moment 
acting on the asteroid due to the sunlight: 

G N 

nd.u) = ^ £/<(«) a) 

a i=i 

c N 

M(d,u) = £5>,x/ 4 (ii) (2) 

a i=i 

where F and M are the total force and moment acting on the asteroid at a 
given instant, ft~lx 10 14 kg km/s 2 is the solar constant, d is the distance 
from the asteroid to the sun, f i is the force acting on each surface facet, both 
at the instant of insolation and after a characteristic time-lag for reemission, 
rj is the position vector from the presumed center of mass of the body to 
the surface facet, and TV is the number of facets for the shape model. When 
expressed as a Fourier series (or as a Spherical Harmonic series in [12]) these 
forces can be represented as: 

r> oo 

F{d,u) = ~J2lMt)™s{\-6)+B l {8)sm{\-6)} (3) 

d 1=0 
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M(d,u) = ^]T[C^)cos(A-#) + £>,($) sin(A-0)] (4) 
d z=o 

where 5 is the solar latitude, computed from sin 5 = sin 2 sm(tu + v), A is the 
solar longitude, 9 is the asteroid rotation angle about its rotation pole, i and 
w are the asteroid orbit elements of inclination and argument of periapsis in 
the asteroid-fixed frame with the z-axis taken as the rotation pole, and v is 
the true anomaly. The Fourier coefficients are shown as being functions of the 
latitude, however they can be further expanded in constant coefficients using 
Associated Legendre Functions [12]. The Fourier coefficients for the moment, 
Ci and Di, are directly computed from the coefficients for the force, where 
we note that the total coefficient equals the sum of the coefficients for each 
facet. Thus we have the relationships: 

N 

C { = 5>,xA| (5) 

i=i 

N 

D l = ^nxBi (6) 

8=1 

where A] and B\ denote the Fourier coefficients at order I for the ith face. 
In these expressions we suppress the time-lag, although it can be included 
into this general form. 

As described in [HI [12], the rotational dynamics of an asteroid subject to 
these torques can be averaged over one rotation about its axis and about one 
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revolution around the sun. Doing so analytically removes almost all of the 
moment coefficients from consideration, and isolates the secular contributions 
of YORP to the asteroid spin rate, u, obliquity, e, and right ascension, a: 

<* = Gl C 0g (e) (7) 

e = %= (8) 

u)I z a 2 y/\ — e 2 

(C hx (e) + D hy (e)) cos(<f) lag ) + (pi, x {e) - Ci, y (e)) sin(0 w ) 

& = - ;°;<^ ( 9) 

- (-Di^e) - ^i,y(e)) cos(0 /a9 ) + (Ci )ie (e) + A, w (e)) sin(0/ afl ) 

where 0; as is the thermal lag angle, and a and e are the asteroid orbit's semi- 
major axis and eccentricity. This form of the secular equations implicitly 
assumes that only the reemission terms are present. 

3 Corrections due to density inhomogeneities 

A common assumption in all previous analyses of YORP effects for given 
shape models, including the theory given above, has been that the body has 
a uniform density distribution, and thus that the are measured from the 
center of figure for a given shape. In the following we continue to assume 
that the moments and forces are computed from the center of figure of the 



body, but now we will also allow for the center of mass and the principal 
axis orientations of the body to be different from the homogeneous density 
case. The theory developed in [TT| IT2] has the common assumption that the 
dynamics are evaluated relative to the center of mass of the system and in a 
principal axis frame. We will continue with these assumptions and, instead, 
will compute the corrections to the constants given in the center of figure 
frame and shift them to the center of mass and principal axes. 

To account for the mass distribution we specify a vector offset from the 
center of mass to the center of figure, denoted as R, and a transformation 
dyad from the figure principal axes to the true principal axes, T (see [11] for 
a brief discussion of dyads as used in this theory). Thus, the new moment, 
accounting for the mass offset, is computed as: 



where we introduce the CM notation for center of mass and CF for center of 
figure. Note that F C f = Fgm, and thus that A^ C f = A-i,cm an d Bi,cf = 
Bi^cm, and so we suppress the CM and CF notation for the force coefficients. 
Carrying these changes over to the moment Fourier coefficients we find: 



M CM 



M CF + R x F CF 



(10) 




C UCF + Rx Ax 



(11) 
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Di,cm = D hCF + RxB l (12) 

To map these vectors into the principal axis frame we just pre-multiply with 
T. 

The corrected equations for the constants in Eqns.[7]-[9]are then computed 

as: 

Co, z ,cm = z-T-[Co,cf + RxAo] (13) 
Cx^cm = x-T-[C liCF + RxA 1 ] (14) 
Ci, v ,cm = y-T-[C li cF + R^A 1 \ (15) 
D 1<X , C M = x-T-[D 1 ^F + Ry<B 1 ] (16) 
D hy ,cM = y-T-[D hCF + Rx BJ (17) 

We should recall that all of the above coefficients are actually functions of the 
obliquity. In keeping with previous practice, we compute these coefficients 
at discrete values of obliquity, although it is possible to define a unique series 
for each coefficient in the obliquity, as explicitly shown in |12j . 

If we assume that the degree of density inhomogeneity is relatively small 
we can develop explicit corrections. Then the magnitude of the center of 
figure offset from the center of mass is small relative to the mean radius of 
the body. The implications of this assumption for the re-orientation of the 



body axes can be seen by decomposing the reorientation dyad using axis- 
angle variables: 

T = cos <pU + (1 — cos (j))aa + sin <pa (18) 

where U is the identity dyad, a is the unit vector pointing along the axis 
of rotation for the re-orientation, is the rotation angle about that axis, 
and the notation a denotes the cross-product dyad applied to the vector a. 
Our smallness assumption for the re-orientation means that both the angle 
is small, and that the axis of rotation is close to the original z-axis, or 
that to first order the axis equals a = a x x + a y y + z. This vector has unit 
magnitude as a% and can be ignored as being small. Ignoring higher-order 
combinations of 0, a x and a y yields the first order result 

T ~ U + </>! + ... (19) 

with the main effect being the shift of the x-y axes about the nominal rotation 
pole. 

Applying these smallness assumptions to Eqns. |T3]fT7| yields the explicit 
corrections: 

Co,z,c*m = Co,z,cf + RxAo : y — R y Ao !X (20) 

Cl,x,CM = Cl,x,CF + RyAi )Z — R z A\^ y 
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+0 [Ci,y,CF + RzM,x - RxM,z\ (21) 
Cl,j/,CM = Cl,2/,CF + R z Ai jX — R x A 1>z 

-0 [Ci, x ,cf + R y A 1>z - R z A hy ] (22) 
Di,x,cm = Di,x,cf + R y B\ )Z — R z B ltV 

+0 [£>i iW) cf + R z Bi, x - R x B l)Z ] (23) 

Dl,y,CM = Dl,y,CF + RzB\ )X — R X B\^ Z 

-0 [D!, x ,cf + R y B hz - R z B 1:y ] (24) 

In the above formulae we have ignored second order combinations of and 
R x and R y . 

Finally, let us consider the moment of inertia about the largest principal 
axis, I z cm- Following from standard results this is computed as: 



Lz,CM 



= z-T ■ [Icf + MR -R\-T T -z (25) 

From the first order dyad correction term we note that z ■ T — z, as the first 
order correction terms are orthogonal to z. Thus we arrive at the updated 
moment of inertia, correct to the second order in the re-orientation terms: 

I z ,cm = I z ,cF-M(R 2 x + R 2 y ) (26) 

If we consider the mass offset terms to be small the moment of inertia is 
essentially unchanged. 
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4 Representative Corrections for Itokawa 

By computing the correction terms related to the force coefficients we can 
gauge the potential strength of this effect independent of a measurement 
of the mass offset. We perform the computation for the asteroid Itokawa. 
We will not focus on the obliquity corrections, as they will not be directly 
measurable, but focus on the spin-rate correction term coefficients Aq jX and 

A ,y 

Shown in Fig. [I] is the YORP spin rate torque coefficient C 0>2i cf and in 
Fig.[2]the averaged force coefficients, all as a function of obliquity (the current 
value of Itokawa obliquity is approximately 178°). In all of these plots we 
show results from a range of surface resolution models. This is to emphasize 
the fact that the YORP effect is sensitive to the asteroid shape, and that 
it evidently has not converged even at the 0.5 meter resolution shape model 
|13j . We only compute the emission terms, as these are expected to be the 
dominant contributor. 

For Itokawa's obliquity of 178° we note that Co, z ,cf, A q<x and A 0tV are all 
negative. Thus Itokawa should currently be experiencing a deceleration in its 
spin rate. A positive shift in the center of mass from the center of figure in 
the x axis (or R x < 0) will cause an increase in the acceleration, making the 
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body decelerate less rapidly for a small enough shift, and a positive y axis 
shift in the center of mass (R y < 0) will cause a decrease in the YORP rate 
of the asteroid, making it decelerate faster. A convenient way to measure the 
strength of this effect is to compute the necessary offset to cause the YORP 
spin rate of change to equal zero. Using the 4 meter resolution model we note 
that a shift of 11 meters in the center of mass from the center of figure along 
the +x-axis or a shift of 15 meters along the — y-axis would cause the YORP 
rate to be zero. For the 0.5 meter resolution model we find that the zero 
line corresponds to a +x-axis offset of 19 meters and a — y-axis offset of 30 
meters. This can be graphically seen by superimposing the zero- YORP rate 
location of the center of mass onto the pole projection of Itokawa, centered at 
the center of figure, shown in Fig. [3] The smallest shift in the center of mass 
that can cause a zero YORP effect is 16 meters for the 0.5 meter resolution 
shape model and less than 9 meters for the 4 meter model. 

In j2j [15] and [8] detailed descriptions of the morphology of the Itokawa 
"neck" region are given which find that loose regolith has migrated into 
this region from neighboring locations with higher elevation. The remaining 
surface of the asteroid is found to be rather rocky and evidently devoid of 
a finer regolith covering. Assuming a common grain density for all of the 
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material, we would expect that the neck region would then have a higher 
bulk density than the surrounding regions. This would cause a shift in the 
center of mass from the center of figure towards the neck region, which is 
consistent with an increase in the YORP deceleration rate, giving it a smaller 
magnitude, and potentially even giving it a positive rotation rate. This also 
provides an explanation for why the YORP deceleration of Itokawa has not 
been detected as of yet. 

5 Discussion and Conclusions 

The inclusion of density inhomogeneities in the Itokawa mass distribution, 
with a concentration of density in the neck region of that body, may account 
for its non-detection, to date, and implies that its detection epoch will be 
delayed. An implication of this result is the possibility to constrain the mass 
distribution of Itokawa once its YORP rotational acceleration is measured. 
Then the center of mass of the body can be constrained to lie on lines similar 
to those shown in Fig. [3j as a function of which shape model is used. Given 
such results it would then be possible to investigate different assumed regolith 
depths in the neck region to constrain the observed shift. 

This result is currently only applicable to the asteroid Itokawa, as it is 
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the only asteroid which should be affected by YORP and for which we have 
a precise shape, spin pole and mass estimate. Future missions are currently 
being planned and proposed to asteroids, some of which should be affected 
by YORP. This theory should be of especial interest, as measurements of the 
shape, mass and YORP effect can then be used to shed additional insight 
into an asteroid's mass distribution. 
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Figure Captions 



Figure 1: YORP torque coefficients Co,z,cf for Itokawa at three different 
surface resolutions. 



Figure 2: YORP force coefficients a) A 0/x and b) A 0jV for Itokawa at three 
different surface resolutions. 



Figure 3: Zero deceleration center of mass line projected onto a pole-on view 
of Itokawa. If the center of mass is to the right of this line, Itokawa will be 
accelerating due to YORP. If it is to the left of this line, it will be decelerating. 
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